Energy dependence of pion in-medium effects on tt /7r+ ratio in heavy-ion collisions 
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Within the framework of a thermal model with its parameters fitted to the results from an isospin- 
dependent Boltzmann-Uehling-Uhlenbeck (IBUU) transport model, we study the pion in-medium 
effect on the charged-pion ratio in heavy-ion collisions at various energies. We find that due to 
the cancellation between the effects from pion-nucleon s-wave and p-wave interactions in nuclear 
medium, the 7r~/7r+ ratio generally decreases after including the pion in-medium effect. The effect is 
larger at lower collision energies as a result of narrower pion spectral functions at lower temperatures. 
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The density dependence of nuclear symmetry energy 
Esymip) is important for understanding isospin effects 
in heavy-ion reactions and properties of rare isotopes as 
well as many phenomena in nuclear astrophysics [H-Q. 
Extensive studies using various probes during past many 
years have significantly constrained the values of £'sym 
at subsaturation densities (for a recent review, see, e.g., 
Ref. Q). The behavior of -Esum(p) at suprasaturation 
densities remains, however, uncertain in spite of many 
studies based on the 7r~/7r+ ratio fsl-fl^. K°/K~^ ra- 
tio n/p or t/'^He rati o [TEI . and the ratio of neutron 
and proton elliptic flows [l6[ in heavy-ion collisions as 
well as the properties of neutron stars |17| . In particular, 
contradictory conclusions were obtained on Esymip) from 
the 7r-/7r+ ratio To describe the FOPI data from 

GSI, results by Xiao et al. Q] based on the IBUU trans- 
port model together with an isospin- and momentum- 
dependent interaction favors a supersoft symmetry en- 
ergy, while those of Feng et al. [ll| based on an im- 
proved isospin-dependent quantum molecular dynamics 
model requires a stiff symmetry energy. A more recent 
study by Xie et al. |12| based on an improved isospin- 
dependent Boltzmann-Langevin approach also favors a 
supersoft symmetry energy. 

In all these studies of the symmetry energy effect on 
7r^/7r+ ratio, medium effects on the pion production 
threshold in nucleon-nucleon collisions that may lead 
to a n op posite effect from that of the symmetry en- 
ergy were neglected. Also, pions were treated as 
free particles in nuclear medium. It is known that due 
to its coupling to the nucleon-particle-nucleon-hole and 
A-particle-nucleon-hole (A-hole) excitations, the disper- 
sion relation of a pion in nuclear medium is softened and 
the strength of its spectral function is enhanced at lower 
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energies [l^l . Because of the isospin dependence of pion- 
nucleon interactions, the tt^ couples more strongly with 
neutrons than the ir'^ and thus has an even softer dis- 
persion relation in neutron-rich nuclear matter than 7r+. 
On the other hand, the pion mass in neutron-rich mat- 
ter is larger for 7r~ than for 7r+ as a result of the pion 
s-wave interaction as shown in Ref. (l9| based on the 
chiral perturbation calculation. Within a thermal model 
for heavy- ion collisions [l^ , effects from the pion-nucleon 
p-wave and s-wave interactions were, however, found to 
largely cancel out, resulting in only a slight reduction 
of the 7r~/7r+ ratio compared to that without the pion 
in-medium effects. 

The pion in-medium effect studied in Ref. (20| was for 
Au-f Au collisions at only one energy of 400 AMeV, us- 
ing 2po {po — 0.16 fm~'^ being the normal nuclear matter 
density), T = 43.6 MeV, and (5iikc = 0.135 for the den- 
sity, temperature, and isospin asymmetry of produced 
matter. Since different densities and temperatures are 
reached in heavy-ion collisions at different beam energies, 
the TT~ /tt~^ ratio is expected to be modified differently by 
pion in-medium effects. It is therefore worth to investi- 
gate the beam energy dependence of pion in-medium ef- 
fects to find the most suitable collision energies where the 
tt^/tt^ ratio is more sensitive to the density dependence 
of nuclear symmetry energy than the pion in-medium ef- 
fects. In this paper, we carry out such a study within 
a similar framework as in Ref. [l^l but using thermal 
model parameters that are more accurately determined 
from the IBUU transport model calculations. 

In high-energy heavy-ion collisions, pions and A reso- 
nances are mostly produced from the initial high-density 
(p > Po) phase. Assuming that this high-density mat- 
ter is in thermal equilibrium with baryon density ps = 
Pn + Pp + Pa + Pa + Pa + Pa^^ isospin asymmetry 
<5iikc = [pn ~ Pp + Pa- - PA++ + Pao/3 - Pa+/3)/pb, 
and temperature T . The densities of nucleons, A reso- 
nances, and pions with isospin states to,-, m-r, and mt 
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can then be expressed, respectively, as 
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Pa 
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In the above, U^j!^^ is the momentum-independent nu- 
cleon mean-field potential with an isoscalar as well as an 
isovector part. The former is fitted to the binding en- 
ergy Eo = —16 MeV and incompressibility Kq = 230 
MeV of normal nuclear matter at saturation density 
Pq. The latter is taken to be from either the super- 
soft symmetry energy that vanishes at about 3po and 
is denoted by a; = 1 or a stiff one that increases al- 
most linearly with density and is denoted by a; = in 
Ref. [21|. These symmetry energies correspond to cur- 



rent uncertainties in -Esym(p) at suprasaturation densi 

ties. Ur = E™.,„J(i- ^ 
A mean-field potential with 
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the 



rirp 1 1 TTit i ttIt- ) being the 
Clebsch-Gordan coefficient from the isospin coupling of 
the A resonance with nucleon and pion. /xs and /ig are 
the baryon and charge chemical potentials, respectively. 
P^^(M) and S'™*(a;,p) are the mass distribution of A 
resonance and the pion spectral function including effects 
due to both pion-nucleon s-wave and p-wave interactions, 
and they are determined from a self-consistent calcula- 
tion as described in details in Ref. [2011. 
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FIG. 1: (Color online) Time evolution of central baryon den- 
sity (top panels), baryon multiplicity (middle panels), and 
average baryon density (bottom panels) in the high-density 
(p > po) phase of central Au-|-Au collisions at different beam 
energies. Horizontal lines indicate the average values used in 
the thermal model. 

The and za in Eq. ([T]) are the fugacity parameters of 
pious and A resonances, respectively, and their values at 
a given temperature depend on the pion and A resonance 
numbers relative to their equilibrium values. Because 
of the short chemical equilibration times estimated from 



the pion and A resonance production rates in the ab- 
sence of in-medium effects, we assume that both and 
ZA are 1 in this case. To determine their values when 
in-medium effects are included, we first determine the 
baryon density, isospin asymmetry, and temperature of 
the nuclear medium from the IBUU transport model cal- 
culations with the same nucleon and A resonance mean- 
field potentials as used in the thermal model but with- 
out the pion in-medium effects. Figure [T] displays the 
time evolution of the central baryon density, the baryon 
multiplicity, and the average baryon density in the high- 
density phase of central Au+Au collisions at three beam 
energies of 200 AMeV, 400 AMeV, and 1 AGeV. Higher 
central densities are reached at higher collision energies 
but last a shorter duration than at lower collision ener- 
gies. Also, the baryon number in the high-density region 
changes with time, and the average values are about 200 
for 200 AMeV, 240 for 400 AMeV, and 280 for 1 AGeV 
as given in Table HI They are, however, similar when dif- 
ferent £'sym(p) are used in the IBUU model. 

TABLE I: Baryon number N'j^, pion-like particle multiplicity 
Ntt, their ratio N-^/Ng, average baryon density (ps), tv' /tv^ 
ratio from the IBUU model with free pions and symmetry 
energy a:: = 1 (a; = 0) for central Au+Au collisions at dif- 
ferent beam energies; and temperature T, isospin asymmetry 
(5iikc, baryon and charge chemical potentials /ib and /ig, fu- 
gacity parameters 2^ and 2^ for pions and A resonances, and 
7r~/7r^ ratio with pion in-medium effects in the corresponding 
thermal model. 





200 AMeV 


400 AMeV 


1 AGeV 




200 


240 


280 




0.52 (0.47) 


9.59 (9.25) 


70.8 (70.2) 


N^/N% (xlO-^*) 


2.60 (2.36) 


40.0 (38.5) 


253 (251) 


(Pb>/Po 


1.3 


1.5 


1.8 


free -k' /-k^ 


4.34 (4.05) 


2.74 (2.63) 


1.79 (1.76) 


T (MeV) 


32.6 (32.2) 


53.3 (52.9) 


90.8 (90.5) 


5likc 


0.152 (0.128) 0.163 (0.136) 0.211 (0.185) 


/iB (MeV) 


919.9 (920.6) 894.3 (895.8) 


814.9 (817.0) 


PQ (MeV) 


-23.6 (-22.5) 


-24.8 (-24.5) 


-23.1 (-23.8) 




0.115 (0.116) 0.190 (0.189) 0.392 (0.391) 


2a 


0.422 (0.417) 0.667 (0.634) 0.813 (0.806) 


in-med. -k' /■k'^ 


4.15 (3.85) 


2.70 (2.59) 


1.80 (1.77) 



The multiplicities of pion-like particles at different col- 
lision energies are also given in Table U and they are 
slightly larger for x = 1 than for x = Q. For the tem- 
perature of the high-density phase, it is determined from 
fitting the ratio of the pion-like particle number to the 
baryon number using the thermal model. Since the av- 
erage baryon density of the high-density phase changes 
with time as shown in the bottom panels of Fig. [1] 
we use the mean values of l.Spo for 200 AMeV, 1.5po 
for 400 AMeV, and l.Spo for 1 AGeV for the thermal 
model calculations. As expected, the resulting temper- 
ature and multiplicity of pion-like particles in the high- 
density phase increase with increasing collision energy as 
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shown in Table ID For the isospin asymmetry Sny^c used 
in the thermal model, it is determined by fitting the fi- 
nal 7r'/7r+ ratio from the IBUU calculations shown in 
Table [J It is seen that (Snkc is larger for a; = 1 , consis- 
tent with the argument that a softer symmetry energy 
leads to a more neutron-rich high-density phase and thus 
a larger 7r~/7r+ ratio as 7r~(7r+) is mostly produced from 
neutron-neutron (proton-proton) scatterings. 

Further shown in Table U arc the baryon and charge 
chemical potentials fiB a-nd fiq determined from the 
chemical equilibrium conditions. Since the multiplicity 
of pion-like particles is much smaller than the total nu- 
cleon number in heavy-ion collisions at energies consid- 
ered here, including pion in-medium effects has negligi- 
ble effects on the extracted values for ps, (Jiikc, and T 
as the dynamics of heavy-ion collisions is hardly modi- 
fied. Indeed, the multiplicities of final pions at different 
collision energies are reasonably reproduced by transport 
model calculations with free pions [1, [111, [l^l , and includ- 
ing pion in-medium effects does not affect much the total 
pion yield [12] ■ Therefore, same values oi fiB and fiq are 
obtained with and without the pion in-medium effects. 
However, due to the softening of pion spectral functions 
and the broadening of A mass distributions from pion 
in-medium effects, the multiplicities of pion-like parti- 
cles would increase if they are assumed to be in chemical 
equilibrium with nucleons. As in Ref. [20| . to reproduce 
the pion and A resonance multiplicities from the IBUU 
model with free pions, their fugacity parameters used in 
the thermal model become much smaller than one if pion 
in-medium effects are included as shown in Table ID Their 
values increase, however, with increasing collision energy 
due to a shorter chemical equilibration time. 
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FIG. 2: (Color online) Pion spectral functions at momentum 
p = rriTr with the supersoft symmetry energy x = 1. 

In Fig. [21 we show the pion spectral functions at a typi- 
cal pion momentum p = 771^^, obtained with the supersoft 
symmetry energy x = 1, in the presence of pion-nucleon 
s-wave and p-wave interactions at different temperatures 



corresponding to different collision energies of 200 AMeV, 
400 AMeV, and 1 AGeV. It is seen that the 7r+ spectral 
function peaks at lower energies compared to that of tt", 
resulting in an increase of the tt"*" multiplicity and a de- 
crease of that of 7r~. Also, the widths of pion spectral 
functions become broader with increasing temperature, 
leading to a larger isospin-depcndent pion in-medium ef- 
fects at lower collision energies. 
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FIG. 3: (Color online) Same as Fig. [2] for the A resonance 
mass distributions. 

The temperature dependence of the pion spectral func- 
tions is related to that of A mass distributions, as dis- 
played in Fig. |3] for the supersoft symmetry energy x ~ 1 
as well. Similar to the pion spectral functions, the A 
mass distribution is broader at higher temperatures. At 
lower temperatures or collision energies, there is a larger 
probability for A"'"+ to have lower masses than for A~ 
as shown in the insets. As a result, the production of 
positive pion-like particles is more enhanced than that of 
negative ones. This isospin-depcndent in-medium effect 
is, however, rather small at higher temperatures. Since 
pions in the IBUU transport model arc produced from A 
resonance decay, the mass distribution of A resonances 
near the nucleon and pion threshold is important for de- 
termining the pion yield in heavy-ion collisions near and 
below the threshold energy. 

The final tt^/tt^ ratios after including both pion- 
nucleon s-wave and p-wave interactions at different colli- 
sion energies with different Esym{p) are given in Table |T] 
and compared to their values using free pions in the left 
panel of Fig. 2] It is seen that the pion in-medium ef- 
fects generally reduce the 7r~/7r+ ratio, leading to the 
need of an even softer symmetry energy for reproduc- 
ing their ratio using free pions. In addition, although 
the 7r~/7r+ ratio at lower collision energies is larger and 
more sensitive to the symmetry energy in the absence of 
medium modifications of the production threshold, the 
pion in-medium effects are also larger, especially at ener- 
gies below the pion production threshold. This is consis- 
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tent with the stronger isospin-depcndent pion in-medium 
effects at lower temperatures shown in Figs. [2] and [31 
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FIG. 4: (Color online) Collision energy dependence of (a) 
7r~/7r^ ratios with and without pion in- medium effects and 
(b) relative 7r~/7r+ ratios from pion-nucleon s-wave interac- 
tion, p-wave interaction, both p- and s-wave interactions, and 
X = with respect to that from a; = 1 in free space from ther- 
mal model for central Au-f Au collisions at different energies. 

To understand the relative effect due to the pion- 
nucleon s-wave and p-wave interactions, we compare in 
the right panel of Fig. |4] the 7r~/7r+ ratios due to the 
symmetry energy effect with those due to only the s-wave 
interaction, only the p-wave interaction, and both inter- 
actions. With only the s-wave interaction, the tt^/tt^ 
ratio is significantly reduced. At the collision energy of 
200 AMeV, this reduction is about 14% and is larger 
than the symmetry energy effect of about 7%. Including 
only the p-wave interaction, the tt~ /tt^ ratio is, on the 
other hand, increased appreciably, reaching about 15% at 



the collision energy of 200 AMeV. After including both 
s-wave and p-wave interactions, the effect of the s-wave 
interaction turns out to dominate over that of the p-wave 
interaction, leading to a 4% reduction in the tt~ /tt^ ratio. 
This value is still comparable to that due to the symme- 
try energy effect. At higher collision energies, the pion 
in- medium effect on the 7r~/7r+ ratio becomes, however, 
very small. 

To summarize, we have found via a thermal model that 
the pion in- medium effects reduce the 7r~/7r+ ratio in 
high-energy heavy-ion collisions compared to that using 
free pious in spite of the cancellation between the pion- 
nucleon s-wave interaction, which reduces the Tr~/'K~^ 
ratio, and the p-wave interaction, which increases the 
7r~/7r+ ratio. Although at lower energies the charged- 
pion ratio is more sensitive to the symmetry energy, 
the pion in-medium effect is also larger, especially at 
collision energies below the pion production threshold. 
Our results thus indicate that to understand quantita- 
tively the symmetry energy effect on pion production in 
heavy-ion collisions, it is important to include the isospin- 
depcndent pion in-medium effects, although this is highly 
nontrivial in the transport model. 
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